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ABSTRACT 


To  address  the  challenge  of  small  event  monitoring  for  the  Comprehensive  Nuclear- Test-Ban  Treaty  (CTBT), 
there  is  great  interest  to  extend  the  successful  teleseismic  mb-Ms  discriminant  to  regional-distance  applica¬ 
tions.  Among  the  outstanding  issues  are  (1)  how  well  can  we  expect  the  discriminant  to  perform  for  regional 
data  as  a  function  of  source  size  and  (2)  what  is  the  physical  basis  behind  its  performance.  In  a  recent  paper 
Patton  (2000),  Mw:  mb  scaling  relationships  for  Pn  and  Lg  waves  with  1-s  periods  were  developed  for  earth¬ 
quakes  and  explosions  located  in  distinct  tectonic  regions  and  geologic  materials  around  the  world.  Among  a 
number  of  findings  of  that  study  is  the  result  that  mb(Pn)  for  explosions  scales  at  a  significantly  higher  rate 
than  mb(Lg)  in  the  Mw  range  ~3.5  -  6.0.  On  the  other  hand,  Pn  and  Lg  scaling  rates  do  not  differ  for  earth¬ 
quakes.  On  plots  of  Ms  versus  mb(Pn),  the  scaling  results  suggest  that  earthquake  and  explosion  populations 
converge  at  Ms  ~  1  -2,  while  populations  may  or  may  not  converge  on  plots  of  Ms  versus  mb(Lg),  and  show 
better  separation  at  small  magnitudes.  Observations  of  regional  Ms  for  small-yield  explosions  at  the  Nevada 
Test  Site  (NTS)  confirm  the  scaling  predictions.  Thus,  mb  based  on  regional  shear  phases  may  serve  as  a  bet¬ 
ter  mb-Ms  discriminant  than  mb  based  on  regional  compressional  phases,  a  surprising  conclusion.  These  find¬ 
ings  certainly  cannot  be  explained  by  simple  theories  of  the  explosion  source,  and  I  propose  a  physical  model 
combining  well-established  scaling  models  and  the  latest  understanding  of  Lg  generation  by  underground 
nuclear  explosions. 
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OBJECTIVE 

The  objective  of  this  project  is  (1)  to  characterize  mb-Ms  scaling  for  earthquakes  and  explosions  using 
regional-distance  observations  and  (2)  to  provide  a  physical  basis  for  the  scaling  and  for  discrimination  perfor¬ 
mance  at  small  magnitudes.  One  of  the  greatest  challenges  of  monitoring  a  CTBT  is  discriminating  man¬ 
made  seismic  events  from  natural  events  at  low  source  strengths,  and  this  project  seeks  an  understanding  of  the 
potential  at  small  magnitudes  using  extensions  of  the  mb-Ms  discriminant  to  regional  data. 


RESEARCH  ACCOMPLISHED 


A  recent  study  of  regional  scaling  relations  found  that  mb(Pn)  scales  at  a  higher  rate  than  mb(Lg)  on  plots  of 
Mw  versus  mb  for  explosions  detonated  at  NTS  (Patton,  2000).  Scaling  relations  for  explosions  at  Semipalat- 
insk  and  other  hard-rock  test  sites  suggested  similar  differences  between  Pn  and  Lg  scaling  rates,  but  the  data 
sets  were  smaller  and  showed  more  scatter  than  NTS  observations,  and  as  a  result,  estimates  of  scaling  rates 
were  less  certain.  In  contrast  to  explosions,  earthquakes  did  not  show  differences  in  scaling  rates  for  Pn  and 
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Lg  waves.  In  this  paper,  I  provide  further  evidence  that  source  scaling  of  1-Hz  Pn  and  Lg  waves  is  different 
for  NTS  explosions,  and  I  propose  a  physical  model  to  explain  why  this  might  be.  The  verification  signifi¬ 
cance  of  this  work  relates  to  characterization  of  mb-Ms  scaling  and  to  the  discrimination  potential  at  small 
magnitudes.  This  characterization  is  the  subject  of  current  research,  but  some  tentative  results  will  be  exam¬ 
ined. 

Both  mb(Pn)  and  mb(Lg)  were  measured  for  most  explosions  conducted  at  NTS  in  order  to  develop  better 
technologies  for  estimating  seismic  yield  (e.g.,  Nuttli,  1986;  Patton,  1988;  Vergino  and  Mensing,  1990).  Of 
special  interest  to  this  project  is  the  frequently  overlooked  fact  that  yield-scaling  reported  for  Lg  and  Pn  waves 
differed  for  explosions  detonated  in  tuff,  as  summarized  in  the  table  below. 

Table  1:  Yield-scaling  results  for  NTS  explosions  (after  Patton;  Vergino  and  Mensing) 


Phase 

Test  Region 

Intercept 

Yield 

Slope 

Gas  Porosity 
Slope 

Lg 

Yucca,  AWT 

3.52 

0.96* 

- 

Lg 

Yucca,  AWT 

4.13 

0.75 

-0.025 

Lg 

Yucca  &  Pahute,  BWT 

4.18 

0.76* 

- 

Pn 

So.  Yucca,  AWT  &  BWT 

3.84 

0.91 

-0.027 

Pn 

Pahute,  AWT  &  BWT 

3.87 

0.91 

-0.027 

*  Yield  slope  when  gas  porosity  is  left  out  of  the  regression  model. 


Patton  (1988)  found  that  after  accounting  for  the  effects  of  gas  porosity,  yield  scaling  of  1-Hz  Lg  waves  for 
explosions  above  the  water  table  (AWT)  is  the  same  as  it  is  for  explosions  below  the  water  table  (BWT).  The 
same  effect  of  gas  porosity  is  seen  for  Pn  waves  (Vergino  and  Mensing),  but  the  slope  on  log  yield,  0.91,  is  sig¬ 
nificantly  greater  than  it  is  for  Lg  waves,  0.76.  Thus,  Pn  scaling  rates  with  official  yield  are  20%  higher  than 
they  are  for  Lg  waves.  For  comparison,  estimates  of  scaling  slopes  on  plots  of  Mw  versus  mb  for  NTS  explo¬ 
sions  are  1.37  ±  0.03  and  1.22  ±  0.03  for  Pn  and  Lg  waves,  respectively  (Patton,  2000),  a  12%  difference. 


Direct  comparisons  of  mb(Pn)  and  mb(Lg)  should  also 
reveal  differences  in  scaling  for  NTS  explosions.  Figure  1 
shows  such  a  plot  for  243  explosions  ranging  ~3  orders  of 
magnitude  in  yield.  A  linear  regression  yields  a  slope  of 
1.13  ±  0.02.  In  other  words,  mb(Pn)  scales  at  a  rate  13% 
higher  than  Lg  waves.  This  estimate  is  about  the  same 
percentage  as  the  one  obtained  from  Mw  versus  mb  data, 
which  used  explosions  detonated  in  a  variety  of  media, 
just  as  Figure  1  does.  The  yield-scaling  results  in  Table  1 
pertain  just  to  explosions  fired  in  tuff  and  rhyolite.  In  any 
case,  summarizing  for  NTS  explosions,  scaling  rates  are 
not  the  same  for  regional  phases  based  on  compressional 
(Pn)  and  shear  wave  (Lg)  phases,  as  various  scaling  rela¬ 
tionships  reveal  that  rates  for  Pn  are  ~13  to  20%  higher 
than  they  are  for  Lg. 


Figure  1.  Regional  mb  for  NTS  explosions.  Dash  line  is 
a  regression  fit  to  the  data;  solid  line  represents  equality. 


A  Physical  Model.  In  this  sub-section,  a  model  combining  the  basic  tenets  of  the  Mueller-Murphy  model 
(Mueller  and  Murphy,  1971;  M-M)  with  recent  insights  into  Lg  generation  by  underground  explosions  is  pro¬ 
posed  to  explain  the  differences  in  source  scaling  between  Pn  and  Lg  waves.  In  this  development,  I  use  the 
concept  of  an  “effective”  source  function.  An  effective  source  function  combines  the  features  of  a  traditional 
source  function  with  the  effects  of  complex  wave  propagation  on  seismic  wave  generation.  In  this  case,  the 
propagation  effects  occur  so  close  to  the  source  that  it  is  impossible  to  separate  them  in  far-held  data  from 
intrinsic  source  properties.  An  example  of  such  propagation  effects  is  near-held  Rg-to-S  scattering.  Gupta  el 
al.  (1992)  proposed  that  this  scattering  mechanism  is  particularly  effective  for  generating  low-frequency 
crustal  S  waves  (<  2  Hz )  by  explosions  because  of  their  shallow  depth  of  burial  and  efficient  Rg  excitation 
characteristics.  Subsequent  research,  both  observational  (Patton  and  Taylor,  1995;  Gupta  el  al.,  1997;  Myers 
el  al.,  1999)  and  synthetic  (e.g.,  Jih,  1993;  1995),  has  supported  this  proposed  mechanism.  The  model  of  Rg- 
to-S  scattering  has  undergone  some  refinements  with  subsequent  research.  In  particular,  excitation  nulls 
observed  in  spectral  amplitudes  of  Lg  waves  (Patton  and  Taylor,  1995;  Gupta  el  al.,  1997)  were  interpreted  to 
be  caused  by  the  imprint  of  Rg  source  spectra  excited  by  a  compensated  linear  vector  dipole  (CLVD).  These 
nulls  are  related  to  focal  depth  of  the  CLVD  source,  which  is  a  kinematic  description  of  spall  phenomena 
occurring  in  the  geologic  strata  over  the  explosion  point  source.  Thus,  while  the  explosion  source  is  an  effi¬ 
cient  source  of  Rg  waves,  spallation,  which  reaches  to  the  free  surface  and  has  a  centroid  depth  even  shallower 
than  the  explosion,  appears  to  be  even  more  efficient  generating  Rg  waves  that  scatter  into  1-Hz  Lg  observed 
in  the  far-held.  Another  key  observation  is  that  effects  of  Rg-to-P  scattering  are  absent  in  the  spectra  of  Pg 
waves  from  underground  explosions  (Patton  and  Taylor,  1995),  suggesting  that  this  mechanism  plays  a  lesser 
role  in  the  excitation  of  regional  P  waves. 

Based  on  this  evidence,  I  argue  that  differences  in  scaling  rates  of  Pn  and  Lg  waves  might  be  caused  by  the 
fact  that  effective  source  functions  for  these  phases  are  not  similar,  nor  do  they  scale  alike.  The  model  I  will 
develop  is  illustrated  in  the  figure  below. 


Log  Frequency  ,  Hz 

Figure  2.  Effective  Pn  and  Lg  source  functions  for  explosions  with  yields  of  100  and  1  kt  (solid  lines)  and 
idealized  Rg-to-S  scattering  functions  (dash  lines).  For  illustration  purposes,  the  Lg  source  function  at  100 
kt  is  simply  the  M-M  spectrum  plotted  for  Pn  waves,  while  for  1  kt,  it  is  a  hybrid  formed  by  vertically  shift¬ 
ing  the  100-kt  spectrum  down,  interpolating,  and  suturing  the  1  -kt  M-M  spectrum  at  ~3  Hz. 


The  Pn  source  functions  are  represented  by  M-M  spectra  for  100  and  1  kt  explosions  in  granite  with  P-wave 
speeds,  a,  of  5  km/s,  S-wave  speeds,  (3,  of  3  km/s,  and  density,  p,  2.5  g/cm3.  As  source  size  decreases,  the 
spectral  overshoot  moves  out  of  the  1-Hz  band  to  higher  frequencies.  Lg  effective  source  functions  are  mod¬ 
eled  as  a  composite  of  the  M-M  spectrum  and  a  narrow  band  spectrum  representing  the  contribution  of  Rg-to- 
S  scattering.  While  the  M-M  spectrum  scales  with  source  size,  the  scattering  spectrum  does  not  because  it  is  a 
property  of  the  near-source  velocity  medium.  The  bump  in  Lg  spectra  is  hypothesized  to  persist  at  small 
yields  for  frequencies  near  1  Hz,  because  Rg-to-S  scattering  is  strong  in  this  band  for  large  and  small  events. 

To  relate  this  model  to  observations  of  scaling  differences,  an  analytical  expression  is  needed  for  the  slope, 
Amb(Pn)/AMw.  First  note  that 

Mw  ~  2/3  •  log  M0  ~  2  •  (log  +  log  Ia)  /  3  , 

and 

AMW  =  2  •  (Alog  'P00  +  Alog  Ia)  /  3  and  Ia  =  pa2  , 

where  M0  is  the  seismic  moment  of  an  explosion,  which  equals  47tIavF00  (Aki  el  al.,  1974),  and  'Poo  is  the  static 
level  of  the  reduced  displacement  potential  (RDP).  Ia  is  referred  to  as  source  medium  “impedance.”  Second, 
using  equation  4.88  of  Aki  and  Richards  (1980),  we  can  relate  the  far-held  amplitude,  A(f),  to  the  source  func¬ 
tion,  S(f),  as  follows 


mb(Pn)  ~  log  A(1  Hz)  ~  log  S(l)  -  1/2  •  log  otIa  , 

and 

Amb(Pn)  =  Alog  S(l)  -  1/2  •  Alog  ala  , 


where  Pn  propagates  as  a  body  wave  in  the  mantle  and  S(l)  refers  to  the  source  amplitude  in  the  1-Hz  band. 
Therefore,  an  analytical  expression  for  the  scaling  rate  of  Pn  waves  is 


Amb(Pn)/AMw  =  3/2  • 


Alog  S(l)  -  1/2-Alog  ala 
Alog  'Poo  +  Alog  Ia 


(1) 


Invoking  the  Mueller  and  Murphy  source  model  and  letting  A  be  lPM  for  a  1  -kt  explosion,  then  AW0'76  is 
for  a  W-kt  explosion.  If  B  is  amplitude  of  the  source  function  in  the  1-Hz  band  for  a  1  -kt  explosion,  then 
BKW0'76  is  the  amplitude  for  a  W-kt  explosion,  where  K  is  a  measure  of  overshoot  in  the  source  function. 
Thus,  in  the  range  of  1  to  W  kt,  the  scaling  rate  can  be  written  in  terms  of  overshoot  parameter  K  and  yield  W, 


Amb(Pn)/AMw  =  3/2  • 


log  K  +  0.76-log  W  -  1/2 -Alog  al0 


0.76-log  W  +  Alog  Ia 

where  Alog  Ia  is  the  logarithm  of  impedance  ratio,  I(XW/I(X1,  and  superscripts  refer  to  explosion  yields. 


(2) 


A  similar  development  is  followed  to  derive  an  expression  for  Amb(Lg)/AMw,  except  the  source  is  a  composite 
of  the  explosion  and  Rg-to-S  scattering.  Equation  (1)  above  is  assumed  to  apply  to  Lg  waves  also,  and  the 
scaling  of  VP00  obeys  the  M-M  model,  just  as  it  does  for  Pn  waves.  The  development  departs  from  that  of  Pn 
waves  for  the  1-Hz  band,  where  B'  is  amplitude  of  the  explosion  source  function,  and  <I>  is  amplitude  from  Rg- 
to-S  scattering.  Patton  (1990)  found  that  the  strength  of  the  spall  source  (which  is  proportional  to  the  mass  of 
spalled  material)  scales  as  W  ,  very  nearly  the  same  scaling  as  'Poo  of  the  explosion.  It  follows  then  that 
amplitudes  for  the  composite  source  are  B'  +  <f>  for  a  1  -kt  explosion  and  [KB'  +  <J)]W°'76  for  a  W-kt  explosion, 
assuming  linear  superposition  of  the  explosion  and  Rg-to-S  sources.  The  scaling  rate  for  Lg  waves  can  be 
written  in  terms  of  W,  K  and  a  new  parameter,  K', 

log  (K+K'-D/K'  +  0.76-log  W  -  1/2-Alog  ala 

Amb(Lg)/AMw  =  3/2  • 


0.76-log  W  +  Alog  Ia 


(3) 


where  K'  =  1  +  <J>  /  B'.  For  K'  >  2,  the  source  amplitude  due  to  scattering  is  as  large  as  or  larger  than  the  ampli¬ 
tude  of  the  explosion  source  function.  In  the  following  sub-section,  I  calculate  the  range  of  acceptable  models 
as  a  function  of  the  material  parameter  Ia  under  the  constraints  imposed  by  measured  scaling  rates  for  Pn  and 
Lg  waves  and  associated  uncertainties. 


Solution  Subspace.  Figure  3  shows  solutions 
for  K  and  K'  in  the  subspace  where  the  cen¬ 
tral  value  of  Am^Pnl/AM^  1 .37,  is  assumed. 
Acceptable  solutions  are  bounded  by  2-0 
errors  on  the  estimate  of  Lg  scaling  slope. 
The  plot  spans  reasonable  values  for  the 
overshoot  parameter  K,  based  on  realizations 
of  M-M  source  functions  for  media  at  NTS 
and  yields  of  -100  kt.  Notice  that  many 
solutions  give  K'  >  2,  indicating  Lg  source 
amplitudes  from  Rg-to-S  scattering  are  larger 
than  source  amplitudes  from  the  explosion. 
However,  many  solutions  have  unrealistically 
large  K'  values,  especially  for  smaller  imped¬ 
ance  ratios,  where  solution  space  shows  rapid 
increases  of  K'  values.  In  order  to  constrain 
solution  space  further,  knowledge  of  an 
appropriate  range  for  the  source  impedance 
ratio  is  needed  for  NTS.  This  is  readily  avail¬ 
able  from  a  large  measurement  data  base  that 
was  developed  and  maintained  by  weapon 
laboratories  for  nuclear  test  containment 
studies. 


K 
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Figure  3.  Acceptable  solutions  (or  source  models)  are  located  in 
the  shaded  region.  The  central  value  on  Lg  scaling  rate,  1.22,  is  the 
locus  of  points  shown  by  plus  symbols. 


Rather  than  using  in-situ  measurements  of  working  point  densities  and  P-wave  speeds,  I  chose  to  compute 
source  impedance  with  a  method  that  gives  a  better  measure  of  what  seismic  waves  “feel”  in  the  source  region. 
Details  were  presented  in  a  poster  at  the  22nd  Annual  Seismic  Research  Symposium  12-15  September  2000, 
but  the  basic  approach  is  to  use  measurements  of  seismic  moment  M0  and  cavity  radii  (from  the  containment 
data  base),  and  calculate  “far-held”  estimates  of  Ia  using  Aki  et  a/’s  famous  equation,  M0  =  47tIavF00.  The 
results  show  relationships  for  Ia  between  explosions  conducted  in  different  media  that  are  expected  (e.g.  esti¬ 
mates  of  r„  for  dry  alluvium  are  much  smaller  than  for  wet  tuff).  An  estimate  of  I,^1  was  found  by  taking  an 
average  of  all  data  for  W  <  5  kt,  and  an  estimate  of  Iaw  was  found  by  taking  an  average  of  all  data  for  50  <  W 
<  200  kt.  The  log  ratio,  Iotw/Iot1 ,  is  0.34,  which  is  plotted  as  a  vertical  dash  line  in  Figure  3.  This  result  is 
probably  a  reasonable  estimate  for  Mw-mb  observations  used  to  determine  scaling  rates  of  Pn  and  Lg  waves. 

For  log  impedance  ratios  near  0.35,  solution  space  does  not  show  large  variations  in  K'  as  the  solution  curves 
are  generally  hat.  Values  of  K'  are  quite  reasonable,  and  the  majority  of  solutions  have  K'  >  2.  This  is  consis¬ 
tent  with  the  hypothesis  that  Rg-to-S  scattering  is  an  important  mechanism  for  generating  1-Hz  Lg  waves  from 
underground  nuclear  explosions.  Under  this  hypothesis,  differences  in  Pn  and  Lg  scaling  rates  are  caused  by 
effective  source  functions  with  dis-similar  shapes  and  scaling  behavior.  The  composite  source  gives  a  shal¬ 
lower  slope  because  the  bump  in  Lg  spectra  from  Rg-to-S  scattering  persists  for  small  yields,  while  it  scales 
out  of  the  1-Hz  band  with  decreasing  yield  in  the  M-M  model  used  for  Pn  waves. 

Implications  for  mh-Ms  Discrimination.  The  teleseismic  mb-Ms  discriminant  is  very  effective,  and  its  physical 
basis  is  well  understood  (e.g.  Stevens  and  Day,  1985).  There  is  much  interest  to  extend  the  discriminant  to 
regional-distance  applications  for  monitoring  smaller  events  than  teleseismic  data  can  record.  The  physical 


model  proposed  in  this  paper  builds  on  our  understanding  of  the  explosion  source  and  regional  phase  excita¬ 
tion.  It  provides  a  basis  for  predicting  regional  mb  scaling  at  small  magnitudes.  Based  on  the  results  for  NTS, 
there  is  reason  to  believe  that,  unlike  the  teleseismic  mb-Ms,  the  best  regional  discriminant  might  be  one  which 
utilizes  amplitudes  of  shear-wave  phases,  such  as  Lg,  for  the  mb,  not  compressional  phases. 

This  project  is  in  the  early  stages  developing  data  sets  to  test  mb-Ms  scaling  predictions  at  small  magnitudes. 
Part  of  this  work  involves  developing  technologies  and  calibrating  regions  for  regional  Ms  determinations; 
efforts  at  universities  and  the  lab  are  well  underway  for  western  China  and  surrounding  areas.  Plots  of 
regional  mb-Ms  for  Pn  and  Lg  waves  are  shown  below  (Figure  4)  as  examples  of  some  preliminary  data  sets  I 
have  compiled. 


Figure  4.  Plots  of  regional  mb  versus  Ms  for  areas  in  western  U.  S.  and  Eurasia.  Pn  data  are  just  from  the  western 
U.  S.,  while  transportability  of  Nuttli’s  mb(Lg)  enables  data  from  both  areas  to  be  plotted  on  the  same  graph  (Pat¬ 
ton,  2000).  Measurements  of  Ms  include  regional  and  teleseismic  data,  and  are  reduced  to  a  common  baseline. 
Solid  lines  for  earthquake  and  explosion  populations  are  predictions  from  Mw  versus  mb  scaling  models. 


I  have  exploited  the  intrinsic  transportability  of  Nuttli’s  mb(Lg)  to  plot  explosion  data  from  NTS,  Semipalat- 
insk,  and  Novaya  Zemlya  and  earthquake  data  from  tectonic  regions  of  North  America  and  southern  Eurasia 
on  a  single  graph.  This  is  not  possible  with  mb(Pn)  because  it  has  been  tied  to  teleseismic  mb(P)  on  a  region 
by  region  basis,  and  it  is  well  known  that  mb(P)  shows  significant  variations  worldwide  due  to  regional  bias. 
The  results  to  date  are  very  preliminary,  but  there  is  some  indication  that  source  populations  are  better  sepa¬ 
rated  at  small  magnitudes  for  Lg  than  they  are  for  Pn,  consistent  with  the  scaling  models. 

CONCLUSIONS  AND  RECOMMENDATIONS 


A  major  conclusion  to  date  is  that  mb(Pn)  and  mb(Lg)  based  on  1-Hz  amplitudes  show  significantly  different 
scaling  rates  for  NTS  explosions.  This  is  supported  by  the  analysis  of  large  numbers  of  regional  mb  observa¬ 
tions  with  respect  to  (1)  yield  scaling,  (2)  M0  scaling,  and  (3)  relative  scaling  of  the  respective  magnitudes 
over  a  wide  range  of  source  size.  A  physical  model  which  explains  different  scaling  rates  draws  on  the 
hypothesis  that  effective  source  functions  for  Lg  waves  are  a  composite  of  two  sources,  the  explosion  and 
near-source  Rg-to-S  scattering,  while  Pn  waves  are  generated  only  by  the  explosion.  Consequently,  the  shape 
and  scaling  behavior  of  effective  source  functions  for  Pn  and  Lg  waves  are  dis-similar.  In  light  of  the  evidence 
that  earthquakes  do  not  show  differences  in  scaling  rates,  I  argue  that  regional  mb-Ms  discrimination  may  per¬ 
form  better  for  mb  based  on  shear-waves  (Lg)  in  contrast  to  traditional  compressional-wave  based  mb. 

A  key  question  is  whether  or  not  P-wave  scaling  differs  from  Lg  at  other  test  sites  around  the  world,  especially 


those  with  hard  rock.  At  this  point,  it  is  difficult  to  say  because  the  existing  data  are  limited  in  quantity  and 
range  of  source  size.  This  question  is  important  to  resolve  because  it  will  address  issues  about  the  variability 
of  regional  mb-Ms  performance  and  further  test  our  understanding  of  its  physical  basis.  The  latest  evidence  I 
have  on  this  score  is  presented  in  the  poster  at  the  22nd  Annual  Seismic  Research  Symposium. 

This  is  report  no.  LA-UR-00-3437  of  the  Los  Alamos  National  Laboratory,  Los  Alamos,  NM  87545. 
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